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In this paper we show, through the study of concrete examples, that, depending on the cosmic 
dynamics of the energy density of matter degrees of freedom living in the brane, Randall-Sundrum 
(RS) brane effects can be important not only at short distances (UV regime), but also at large 
cosmological scales (IR regime). Our first example relies on the study, by means of the dynamical 
systems tools, of a toy model based in a non-linear electrodynamics (NLED) Lagrangian. Then 
we show that other, less elaborated models, such as the inclusion of a scalar phantom field, and 
of a tachyon phantom field - trapped in the brane - produce similar results. The form of the RS 
correction seems to convert what would have been future attractors in general relativity into saddle 
points. The above "mixing of scales" effect is distinctive only of theories that modify the right- 
hand-side (matter part) of the Friedmann equation, so that, for instance, DGP-brane models do not 
show this feature. 

PACS numbers: 04.20.-q, 98.80.-k, 98.62.En, 98.80.Cq, 98.80.Jk 



I. INTRODUCTION 



Braneworld models of the Universe, where standard 
model (SM) particles are trapped in a three-space, while 
gravity might propagate also in the extra-space, have 
been the focus of attention in the current decade in the 
search for answers to several outstanding problems in the- 
oretical physics and in cosmology. The Dvali-Gabadadze- 
Porrati (DGP) brane model for instance, was ex- 
plored in an attempt to understand the mystery under- 
laying the present accelerating pace of the cosmic ex- 
pansion, without invoking the mysterious "dark energy". 
The DGP brane model modifies general relativity (GR) 
at large cosmological scales or, in other words, it car- 
ries infra-red (IR) modifications of GR. In this model 
the present speedup of the cosmic expansion is explained 
as an effect of the extra-dimension: at large cosmolog- 
ical scales gravity leaks into the extra space, resulting 
in a weakening of the gravitational bounds of the cos- 
mic structure.^ Acting on the opposite end of the en- 
ergy/curvature spectrum is the Randall-Sundrum (RS) 
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model [1, 0], which carries ultra-violet (UV) modifica- 
tions of GR. This time, however, the motivation was to 
seek for an alternative explanation to the hierarchy prob- 
lem (RSI model 0), and to look for an alternative to 
the Kaluza-Klein compactification procedure (RS2 model 
Q)- The model known as RS2 found important applica- 
tions in cosmology [Sj, in particular for the inflationary 
scenario [6]. In this model a single codimension-1 brane 
with positive tension, is embedded in a five-dimensional 
anti-de Sitter (AdS) bulk spacetime, which is infinite in 
the direction perpendicular to the brane. As said, the 
standard model matter degrees of freedom are confined 
to the brane, while gravitation can propagate in the bulk. 
However, in the low-energy/IR limit, due to the curva- 
ture of the bulk, the graviton is confined to the brane, and 
standard (four-dimensional) GR laws are recovered. The 
latter result can be better understood if we invoke a cos- 
mological application j^: let us consider a homogeneous 
and isotropic Friedmann-Robertson- Walker (FRW) met- 
ric with fiat spatial sections. The Friedmann equation 
on the RS2 brane then reads (we neglect the " dark radi- 
ation" and cosmological constant terms): 

3ij2 = p(l + p/2A), 

where p is the energy density of the matter degrees of 
freedom trapped in the brane, and A is the brane ten- 
sion. The latter is fixed by the curvature scale of the 
bulk, which is constrained to be less than 0.1 mm by 
current tests of Newton's law 0- At very high energy 
density (p A), due to the brane effects, the above 
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equation is fundamentally modified: H on p. As the cos- 
mic expansion proceeds the energy content in the brane 
dilutes and, eventually, as long as p ^ A, the standard 
GR behavior H"^ oc p is recovered. 

It is quite obvious that the latter result is highly de- 
pendent on the cosmic dynamics of the energy density p 
itself. Actually, for fields whose energy density increases 
in the course of the cosmic expansion, the brane effects 
can be important at late-times/large cosmological scales 
also. This simple and obvious fact - when analyzed in 
the light of the cosmological dynamics of fields trapped in 
the brane - seems to contradict the spirit of the Randall- 
Sundrum brane scenario, which is usually assumed to 
modify gravity at high UV energy/curvature, with spe- 
cial impact on primordial (early) inflation scenario. Our 
trivial finding implies that, if the energy density of some 
(perhaps exotic) component of the cosmic mixture does 
not dilute but grows as expansion proceeds, then the cor- 
responding RS-brane model can have implications not 
only for early-time cosmology, but it can have apprecia- 
ble impact also on the fate of the cosmic expansion. 

That is why, in the present paper we will explore the 
impact of RS brane effects on the cosmic dynamics, when 
the energy density of - at least one of the components of ~ 
the matter degrees of freedom trapped in the brane grows 
with the course of the cosmological expansion, instead of 
diluting itself. To this end we will explore the asymp- 
totic properties of concrete models of very different na- 
ture, starting with a very elaborated toy model based on 
a non-linear electrodynamics (NLED) Lagrangian whose 
asymptotic structure is very rich. By investigating how 
the brane effects modify the asymptotic structure of the 
latter model, we will show that RS2 brane effects might 
convert what would have been future attractors in GR 
into saddle critical points. This result will be corrobo- 
rated through the study of other examples of exotic mat- 
ter trapped in the RS brane: i) a scalar phantom field, 
and ii) a tachyon phantom field. The physical conse- 
quences of the above curious ("trivial") finding will be 
discussed. 



II. NLED-BASED TOY MODEL 

As already mentioned, we start by exploring the dy- 
namics of a non-linear electrodynamics (NLED)-based 
toy model. The four-dimensional (4D) Einstein-Hilbert 
action of gravity coupled to NLED is given by: 

S = J d^xy/^[R + L^ + L{F,G)], (1) 

where R is the curvature scalar, L-y ~ the background 
perfect fluid's Lagrangian, and L{F, G) is the gauge- 



invariant electromagnetic Lagrangian, which is a func- 
tion of the electromagnetic invariants F = F^^'^F^i, and 
G = ^eap^uF""^ Q. As usual, the electromagnetic 
tensor is defined as: F^i, = d^A^ — d^A^. Standard (lin- 
ear) Maxwell electrodynamics is given by the Lagrangian 
L{F) = — F/4. The corresponding field equations can be 
derived from the action by performing variations with 
respect to the spacetime metric gfj,u, to obtain: 

where 

T/^u^' = g^AHF) - GLg] - 4F^o.F,"Lf, 

with pj ~ p-yit), p-y — P-yit) - the energy density and 
barotropic pressure of the background fluid, respectively, 
while Lp = dL/dF, Lpp = d'^L/dF^, etc. Variation 
with respect to the components of the electromagnetic 
potential yields to the electromagnetic fleld equations. 

In order to meet the requirements of homogeneous and 
isotropic cosmologies (as, in particular, the one associ- 
ated with FRW spacetimes), the energy density and the 
pressure of the NLED field should be evaluated by aver- 
aging over volume. To do this, we define the volumetric 
spatial average of a quantity X at the time t by (for 
details see [1] and references therein): 

X EE lim — / d^xJ^ X , 

where V = \ d^Xy/—g and Vb is a sufficiently large time- 
dependent three-volume. Following the above averaging 
procedure, for the electromagnetic field to act as a source 
for the FRW model we need to impose that (the Latin 
indexes run over three-space); Ei — 0, Bi = 0, EiBj = 0, 
and also SEiEj ~ —E^gij, SBiBj ~ —B^gij. Addition- 
ally it has to be assumed that the electric and magnetic 
fields, being random fields, have coherent lengths that 
are much shorter than the cosmological horizon scales. 
Under these assumptions the energy-momentum tensor 
of the electromagnetic (EM) field - associated with the 
Lagrangian density L = L(F, G) -, can be written in the 
form of the energy-momentum tensor for a perfect fluid: 

T)f/^ = {PB + Pb) U^Uy - PBgt,u, 

where ps = —L + GLq — ALpE"^, and ps = L — GLq — 
4:{2B^ — E'^)Lp/3, with E and B being the averaged elec- 
tric and magnetic fields, respectively. In what follows, to 
simplify the analysis, we shall consider a FRW universe 
filled with a "magnetic fiuid", i. e., the electric compo- 
nent E will be assumed vanishing. Even this simplified 
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picture can give important physical insights. We will be 
focusing in the study of the asymptotic properties of a 
cosmological model with interesting features, namely a 
phase of current cosmic acceleration and the absence of 
an initial singularity, which was proposed in [s!] (see also 
0) and is based upon the following Lagrangian density: 

L = -^F + aF^ + pF-\ (3) 

where a and /3 are arbitrary (constant) parameters. As 
seen this Lagranagian contains both positive and nega- 
tive powers of F. The second (quadratic) term domi- 
nates during very early epochs of the cosmic dynamics, 
while the Maxwell term (first term above) dominates in 
the radiation era. The last term in ([3]) is responsible for 
the accelerated phase of the cosmic evolution The 
above Lagrangian density yields a unified scenario to de- 
scribe both the acceleration of the universe (for weak 
fields) and the avoidance of the initial singularity, as a 
consequence of its properties in the strong-field regime. 
Recalling that we are considering magnetic universes, 
i.e., F — 2B^, where is an averaged value of the 
magnetic field, ^ the stress-energy tensor associated with 
can be written in the form of an equivalent perfect 
fluid stress-energy tensor with energy density and para- 
metric pressure: ps = 5^(1 — SaB"^ — [3B~'^)/2, and 
PB = B^{1 - 40aB^ + 7l3B-^)/6, respectively. ^ 

In this paper we assume the background fluid to be 
dust cold dark matter (CDM), so that 7 = 1. Our goal 
will be to put the corresponding cosmological equations: 

3iy^ = pcdm + PB, 2H = -pcdm ~ (PB +Pb), 
Pcdrn + SHpcdm = 0, pB + 3H{pB + Pb) = 0, (4) 

in the form of an autonomous system of ordinary differ- 
ential equations (ODE). For this purpose we choose the 
following phase space variables: 

_ PB _ 16aB-^ _ 4:13 

After this choice, the following autonomous system of 
ODE is obtained out of (g]): 



^ For details of the averaging procedure consult Q . 

^ Notice that, for large values of the NLED field, positivity of 
energy requires that B < l/\/8a, while, for small enough values 
of B <SC 1, if one considers positive /3 > 0, then positivity of 
energy implies that B > (7/3)^^*. The existence of the lower 
bound, at first sight might appear problematic, however, given 
that the observational data constraints the parameter /3 to be 
\/\fi\ 4 X 10~^* g cm~^ j9|, then the lower bound on B can 
be admitted without going into conflicts with observations. 



x' = [x- l){x -y + z), y' = -i/(5 -x + y- z), 
z' ^ z{7 + X - y + z). (5) 

Here, for generality of the analysis we shall consider arbi- 
trary a S 3? and /3 G 3?. In what follows in this paper the 
tilde will denote derivative with respect to the new time- 
variable T = In a (c?T = Hdt) , also known as the number 
of e-foldings. The phase space relevant to the present 
study is then given by the following region in (x, y, z): 

^'c; = {(a;,y,^)|0<a;<l, {y,z)e^^, 8x + 2y + z>0}, 

where we have considered the fact that B^/QH^ — x + 
y/A + z/8. In what follows we use the standard deflnition 
of the dimensionless density parameter for the compo- 
nent X: Qx = Px/3i?^, while lux — Px/px will denote 
the effective equation of state (EOS) parameter of the 
component X. Four equilibrium points of ([5]) in are 
found: 

1. Radiation-dominated phase: 

P™d = (x,2/,z) = (1,0,0), nB = i, n^dm^o. 

This is a decelerating expansion solution {q = 1), 
that is fueled by standard radiation with ujb = 1/3. 
The eigenvalues of the linearization matrix corre- 
sponding to this equilibrium point are: Ai = 1, 
A2 ~ —4, A3 — 8, so that it is a saddle in ^j/- 

2. Infra-red NLED-dominated solution: 

n!i^d = (1,0,-8), 17b = 1, ^cdm^O. 

This is a late-time, super-inflationary solution (q — 
—3), where the NLED fluid mimics phantom be- 
havior (ujb = —7/3). This solution exist only for 
z < (negative /3 < 0). In this case there is 
no lower bound on the magnitude of the magnetic 
held. The eigenvalues of the linearization matrix 
for P^Ha are: Ai = -8, A2 = -7, A3 = -12, so 
that this solution corresponds to a late-time (fu- 
ture) attractor. Since, the NLED-magnetic fleld 
mimics phantom behavior, the late-time attractor 
might be associated with a cosmological singularity 
(most probably a big-rip type of singularity) . 

3. Ultra-violet NLED-dominated phase: 

A^rd^ (1,-4,0), nB^i, n,d„.^o. 

This solution corresponds to an early-time, super- 
stiff- fluid solution (wb = 5/3), which is associated 
with super-decelerating expansion (g = 3). This 
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solution exist only for ?; < 0, i. e., if the con- 
stant a is a negative quantity (a < 0). Curi- 
ously, this case does not meet the conditions for 
a bounce (there is no upper bound on the magni- 
tude of the magnetic field), and the corresponding 
cosmology starts with a big-bang singularity, since 
= 0, B ^ 0, ^ H ^ oo. This solu- 
tion is a past attractor in VPc/, since the eigenvalues 
of the corresponding linearization matrix: Ai = 5, 
A2 = 12, A3 = 4, are all positive quantities. This 
means that this point is the starting point of every 
probe path in the phase space 

4. CDM-dominated solution: 

-Pcdm = (0,0,0), CtB^O, n,d,n^l- 

This phase of the cosmic evolution is characterized 
by decelerated expansion {q = 1/2). The EOS pa- 
rameter for the magnetic field is undefined in this 
case. The existence of this solution is necessary for 
the formation of the observed amount of structure. 
This is also a saddle equilibrium point in ^[/, since: 
Ai = -5, A2 = 7, A3 = -1. 

For positive definite a and /3, only the radiation- 
dominated equilibrium point Pnied, and the CDM- 
dominated solution Pcdm, are found in ^[/. The above re- 
sults confirm our expectation that a combination of posi- 
tive and negative powers of the electromagnetic invariant 
F in the NLED Lagrangian, can drive a very interesting 
cosmological scenario leading to accelerated expansion at 
late times. Unfortunately, early-time (primordial) infia- 
tion can not be obtained in the present model. 



The variable v controls the brane regime. 



2{l-v)' 



so that, for instance, the GR-limit (formal limit A — ^ 00) 
corresponds to u = 1. The following autonomous system 
of ODE can be derived out of ([5]): 



x' = Z{l-v)x + 



x~ 1 



{x-y + z), 



y = y 



-5 + (1 - 1;) + (j—^ ] (x-y + z) 

7 +(1-1;)+ ^ ]{x-y + z) 



(1 — v){3v + X — y + z). 



(7) 



At w = 1 the first three equations above coincide with the 
equations ([S]) , which hold for GR with a NLED-magnetic 
field. The phase space of the model can be defined as 
follows: 



brane 



{{x,y,z,v)\ < x < 1, {y,z) e Sft^, 
8x + 2y + z > 0, 0< 1; < 1}. 



The critical points of the autonomous system of ODE 
([7]) in the phase space y^/^jane^ their physical properties 
and stability are discussed below. As in the general rel- 
ativity case, four equilibrium points are found: 



A. NLED on the RS brane 

Now we explore the possible effect of RS2 braneworld 
gravity on the above picture. The corresponding (modi- 
fied) FRW cosmological equations for a RS2 brane with 
CDM and a perfect fiuid of NLED trapped in it, can be 
written as (here we omit the " dark radiation" term) ; 



= Ptot (1 + , Ptot = Pcdm + PB, 



= -{pcdm + PB +Pb) [1 
Pcdm — —3Hpcdm, B = —2HB. 



Ptot \ 

A ) ' 



(6) 



It is convenient to introduce the following phase space 
variables: 

PB 16aB4 4/3 Ptot 



3H^B^ 



1. CDM-dominated solution: 



Prdrr 



(0,0,0,1), ris-o, n,dm = i 



Since q = 1/2, this solution is associated with decel- 
erated expansion. The NLED-magnetic field EOS 
parameter iog is undefined. This is a saddle equilib- 
rium point in vl/^<^"'=. Actually, the eigenvalues of 
the linearization matrix corresponding to this point 
are: Ai = —5, A2 = 7, A3 = —3, A4 — —1. 

2. Radiation-dominated solution: 



P™d - (1,0,0,1), ns^i, ^cdn 







3i?2 



It is also a decelerated-expansion solution {q = 1) 
driven by standard (Maxwell) radiation [ujb ~ 
1/3). As the CDM-dominated phase, this solution 
also represents a saddle equilibrium point in v]/^""'^, 
since the eigenvalues of the linearization matrix are 
of opposite signs: Ai = 8, A2 = 1, A3, 4 = —4. 
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3. UV NLED-dominated solution: 



pUV 
nled 



(1,-4,0,1), Qb^i, n 



cdm 



. 



This solution shares many properties with its sim- 
ilar GR-solution: it is a super-decelerated {q — 3), 
super-stiff state {ujb = 5/3), associated with a big- 
bang- type singularity (H — > oo). The new feature 
here is that the stability properties have been mod- 
ified by the brane effects. Actually, the eigenvalues 
of the linearization matrix in the present case are: 
Ai — 12, A2 = 4, A3 = 5, A4 = —8, so that this 
is a saddle equilibrium point in '^/'^"■'^^ (its similar 
GR-solution is a past attractor). 



4. IR NLED-dominated solution: 



piR 

nled 



(1,0,-8,1), ili 



1, a 



dm 



. 



As in the former case, this solution shares many 
properties with its GR-similar: it is a super- 
accelerated {q — —3), phantom-like solution (a;^ = 
—7/3), possibly associated with a big-rip-type sin- 
gularity. Other properties, the stability in partic- 
ular, have been modified by the brane effects. Ac- 
tually, the eigenvalues of the Jacobian matrix are 
Ai = —12, A2 = —8, A3 = —7, A4 = 4, so that it is 
also a saddle equilibrium point in the phase space 

^Tibrane 

The first thing that is worthy of mention, is the fact 
that all of the critical points found represent saddle equi- 
librium points in the phase space ^p^""^. Note that, 
since v = 1, the four equilibrium points are associated 
with GR. In fact these coincide in almost all aspects with 
the ones found previously for the GR case. The first two 
points Pcdm, and Prad, show no fundamental differences 
with their GR-similar. However, the NLED-dominated 
solutions PniYd ^n/ed' have different stability proper- 
ties than their corresponding GR-solutions: while in the 
latter case P^iYd attractor and P^;^^ was the 

future attractor, in the present case both are saddle equi- 
librium points as already said. This means, in turn, that 
the space-time singularities associated with these critical 
points (big-bang and big-rip singularities respectively), 
might be evaded in the present case. Modification of 
the stability properties of the UV solution was expected 
since, as already mentioned, RS brane effects are appre- 
ciable at high energies/short distances (early times). Ac- 
tually, only at very high energies can the graviton ac- 
quire large momenta along the extra dimension and may 
escape into the bulk (5D) spacetime. The surprise was 
the IR solution: it is expected that at low energies (large 
cosmological scales) the RS2 brane effects can be safely 
ignored. However, while making such statements one has 



to be careful. In the cosmological context, the most ap- 
preciable RS2 brane effect is to modify the Friedmann 
equation: iH^ — ptot{^ + Ptot/'^X)- Hence, at very high 
energy density ptot ^ A (much bigger than the brane 
tension), the Friedmann equation is fundamentally mod- 
ified oc Pjoj. If in the course of the cosmic expansion 
the total energy content of the universe dilutes, then as 
long as Ptot ^ A, one recovers standard GR-Friedmann 
behavior. Now look at the Lagrangian density for the 
NLED-niagnetic field ([3]) considered here. Note that as 
the expansion proceeds the magnetic field F (x di- 
lutes, and the component oc /3_B~^ in ([3]) grows with- 
out limit. This means that the total energy content of 
the universe starts growing at the expense of the NLED 
component so that, at late times, eventually, ptot might 
become much larger than the brane tension once again, 
rendering the brane effects important at late times also. 
It is then to be expected that similar IR modifications will 
be produced if one considers other less elaborated models 
such as, for instance, a scalar phantom field trapped in 
the RS2 braneworld. 



III. THE CASE OF SCALAR/TACHYON 
PHANTOM FIELDS 



Consider now that the source of the Einstein's field 
equations are a background of dark matter and a scalar 
phantom field (f>. In this case the cosmological equations 
coincide with Q but with the replacement of ps, Pb by: 



(8) 



where V{(j)) is the self-interaction potential. To write 
the corresponding cosmological equations in the form of 
an autonomous system of ODE, it suffices to choose the 
following phase space variables: 



The following autonomous system of ODE is obtained: 



V 



(9) 



y 



2 V 
2 V 



— 3x — X 



IP_ 



xy - y' 



(10) 



where H' / H — —3(1 — x^ — y^^jl. As before we can write 
the relevant parameters in terms of the variables y: 
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TABLE I: Critical points of the autonomous system of ODE (|14|l for the potential V{(j>) ~ Vo(j> ^ and their properties. These 
critical points always exist. The parameters Ai and A2 are the eigenvalues of the linearization matrices corresponding to each 

one of the critical points. We have defined: y* — u (A^ + VA* + 36)/6. 



C. Point 


X 


y 




UJ4> 


9 


Ai 


A2 


M 








1 




1/2 


-3 


3/2 


TP 


-Ay*/V3 


y* 


1 


-1 - A2y2/3 


-1 - 




-3 - X^yl 



n, 



1 



y 



y 



Notice, also, that for the deceleration parameter one has 
q ~ —{l + H'/H). The phase space relevant for this case 
is defined by the following set: 



ph 



{{x.y)\0<y^ -x" < 1, y>Q] 



For simplicity we shall consider an exponential self- 
interaction potential for the phantom field: V{(j)) — 
Voexp(— A0), where A is a constant. For this potential 
the system of ODE (fTO|) is closed. Only two equilibrium 
points are found in this case: 



1. Matter-dominated point: 

M = (0,0), n^dra = 1, 



il^ = 0, q = 



2' 



while the equation of state is undefined. The 
eigenvalues of the corresponding linearization (Ja- 
cobian) matrix are: Ai = —3/2, A2 = 3/2, so that 
this a saddle equilibrium point. 

2. Scalar phantom-dominated solution: 



5P = (-A/n/6, + AV6), 



Notice that a;0<— 1, g<— Iso that this is 
a supper-accelerated, phantom-dominated solution, 
independent of whether A is a positive or a nega- 
tive quantity. This critical point is a future attrac- 
tor in the phase space since both eigenvalues 
of the Jacobian matrix are negative in this case: 
Ai=-3-AV2, A2 = -3-A2. 

Lets turn our attention now to the phantom tachyon 
field, for which, instead of the definition ([8|) for the effec- 
tive energy density and pressure, one has: 



V 



(11) 



1 



Notice that the sign of the kinetic energy under the 
square root has been changed. The Klein-Gordon equa- 
tion for the phantom tachyon is 



1 



d4,V 
V 



1 



0. 



(12) 



In order to write the cosmological equations in the form 
of an autonomous system of ODE, in this case it is rec- 
ommendable to work with the following phase space vari- 
ables [13: 



(13) 



V 



Notice that, while the phase space variables © were nor- 
malized by the Hubble parameter H, the ones in (fT3|) are 
not normalized (at least x in (jl3p is not a normalized 
variable). After the above choice, one can obtain the 
following autonomous system of ODE: 



x' = {x^ + 1) 



-3a; + V3' 



y 



where 



y ^y 



H' 



3 d^V 



li 



(14) 



H 



Ul-yWx^ + l). 



In terms of the variables x, y in (jl3p . one has: fl^ = 
y^/Vl + x^, while for the phantom tachyon EOS param- 
eter, W0 = —l — x^. The phase space spanned by the 
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variables x and y, that is relevant for an expanding Uni- 
verse, can be defined as: 



{{x, y) : -1 < a; < 1, < < 1 - x^}. 



In general the system (|14l) is not a closed system of ODE, 
unless one considers the specific case when A = dij)V/V'^^'^ 
is a constant pLO], so that = Vb0~^ (Vo = 4/A^). In 
what follows, for simplicity, we shall consider, precisely, 
the latter kind of self-interaction potential for the phan- 
tom tachyon. The critical points of in this case, as 
well as their relevant properties, are shown in Tab. HI 
As in the former case there are found only two critical 
points: 

1. Matter-dominated solution M — (0,0). This is 
a decelerated expansion-solution and represents a 
saddle point in the phase space ^t- 

2. Tachyon phantom-dominated solution TP = 
(— Aj/,/-\/3, y*)- In this case, since both cj^ < —1, 
and q < —1, it is a super-accelerated, phantom so- 
lution. This equilibrium point is a future attractor 
in '^T- 

The above discussed asymptotic structure is clearly il- 
lustrated in the Fig. (TJ where the phase portraits for 
the scalar phantom field (left-hand panel), and for the 
tachyon phantom field (right-hand panel), respectively, 
are drawn for a given value of the brane tension (A = 2). 
It is apparent how the phase trajectories approach to the 
future attractor: the phantom-dominated solution. 




A. Brane effects 

Lets check now how the brane effects modify the dy- 
namics of the scalar/tachyon phantom field. In this case 
we have to replace p b and p b in the cosmological equa- 
tions (jni) for matter trapped in the brane, by the cor- 
responding expressions for the effective energy density 
and pressure for the scalar phantom field ([5]), and for 
the tachyon phantom field (llip . respectively. In order 
to transform the corresponding field equations into an 
autonomous system of ODE one can use the same vari- 
ables X, y defined in © for the scalar phantom field, or 
in (jl3p for the tachyon phantom field, plus an additional 
variable: 



Ptot 

MP' 



(15) 



In terms of this variable 



Plot 

A 



2(1 -z) 



< z < 1. 



FIG. 1: Phase portraits for the scalar pliantom field (left- 
hand panel), and for the tachyon phantom field (right-hand 
panel) , respectively. The brane tension A has been arbitrarily 
set equal to A = 2. In both cases there are no past attractors, 
while the phantom-dominated solution is always the future 
attractor. 



This means that general relativity is recovered at the 
value z = 1. Notice that at the point z — Q the ratio 
Ptot I ^ is undefined, so that the plane set {(x, y, z)] 2 = 0} 
has to be removed from phase space. However, since we 
are interested in the asymptotic structure, we have to 
study, separately the limiting case z — > 0. 

For a scalar phantom field trapped in the RS brane the 
autonomous system of ODE that can be obtained out of 
the corresponding cosmological equations is the following 
(compare with ([TU)) . V — Vq exp {~X(j))): 
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TABLE II: Critical points of the autonomous systems of ODE (|16() . and (|17|) . and their properties. These critical points always 
exist. The parameters Ai, A2, and A3 are the eigenvalues of the linearization matrices corresponding to each one of the critical 
points. The upper half of the table is for the scalar phantom field, while the lower half is for the tachyon phantom field. We 

have defined: y, = J {X'^ + V>^T36)/6. 



C. Point 


X 


V 


z 








g 


Ai 


A2 


A3 


E 

















undef. 


2 


6 


3 





M 








1 


1 





undef. 


1 
2 


3 
2 


3 
2 


-6 


SP 


A 




1 





1 




\ 2 


-3- A^ 


\ 2 

" 2 


2A2 


E 

















-f 


2 


3 


3 




M 








1 


1 





-1 


1 

2 


-3 


3 

2 


-3 


TP 


Ay. 
V3 


y* 


1 





1 


^ 3 


,2 2 

1 A I/, 
2 




-3 - \\l 


\^yl 



3, H' 



z' = 6(l-z)(z-a;2-y2), 



(16) 



where 



H' 2, f2- z 



H 



{z-x -y ). 



As before (scalar phantom field, no brane), for the phan- 
tom field dimensionless energy density parameter, and 
for the EOS parameter, one has 



ri^ = - x^, 



x'^ + y"^ 
x"^ — y"^ 



Hovifever, the expressions for the CDM dimensionless en- 
ergy density parameter, and for the deceleration parame- 
ter, differ from the ones that appear after equations (fTOt : 



^cdm = z + x^ -y'^, g = -1 + ^ ^ ^ 

The phase space relevant for this case (recall that we 
consider only expanding Universes so that y > 0) is given 
by the following set: 

^'';r' = {{x,y,z)\y^-x^>Q, y > 0, < z < 1}. 

In a similar fashion, for the tachyon phantom field on 
the RS brane we obtain the following autonomous system 
of ODE (compare with V = Vb-^^^): 



x' = ~{x^ + l)(3x- TSAy), 




z' ^ 3(1 - z){z - y2Va?2+T), 



(17) 



where 



H' _ 3 /2-z 



(M) (^-y'^^^)- 



The cosmological parameters of observational interest are 
given by the following expressions: 

2 2 



\Jx^ + 1 ■ 



while, as customary, for the deceleration parameter: q = 
— (1 + H' /H). The phase space where to look for the 
equilibrium points of (jl7|) is 



^brane ^ {{^x,y,z)\ Z\/ x'^ + 1 > y > 0, 0<Z<1}. 



1. Equilibrium points and their properties 

A simple inspection of equations (|16p and ([TTl) , reveals 
that there can be a critical point lying on the phase plane 
(a;,y, 0), which has been deliberately removed from the 
corresponding phase spaces. Actually, the point (0, 0, 0) 
seems to be a critical point of each one of the autonomous 
systems described in both cases. Note, however, that at 
this point the expression for H' / H is undefined, and so 
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-2-10123 
X 







n 





-2-10123 
X 



FIG. 2: Projections of the phase space trajectories for the 
scalar phantom field trapped in the RS brane, onto the phase 
planes (a;, z) and (y, z), respectively. The brane tension A has 
been arbitrarily set equal to A = 2. Notice that not only the 
past asymptotic structure has been radically modified (there 
is now a past attractor: the empty Misner Universe), but also 
the future asymptotic has been modified: there is now no 
future attractor. 



are first and second equations in and second equa- 
tion in ([T7|) . However this fact is a consequence of the 
chosen phase space variables. We can apply a secure ap- 
proach to study the nature of the point (0, 0, 0). In fact, 
we can expand the ODEs in (|16I17|) in the neighborhood 
of (0, 0, 0): a; -I- e^;, 2/ -> + Ej,, z -I- e^, where 
(e^, Cy, Ez) are the small perturbations. If we retain only 
the linear terms in the perturbations, i. e., if neglect 
terms ~ ~ ~ ^x^y ^ ^x^z ~ <^y<^z ^ 0, and 
smaller, then the autonomous systems of ODE ([TS]). and 
([T7)) . transform into the following systems of linear ODE: 




X 




-10 12 
X 



FIG. 3: Projections of the phase space trajectories for the 
tachyon phantom field trapped in the RS brane, onto the 
phase planes (x^z) and {y,z), respectively. The brane ten- 
sion A has been arbitrarily set equal to A = 2. The empty 
Universe (0, 0, 0) is the past attractor for any probe path in 
the phase space, while there is no any future attractor. 

and 

^x ~ N/SAey — Se^;, e'y — 3ej^, = Se^, 

respectively. Hence, in the case of the scalar phantom 
field trapped in the RS brane, the point (0,0,0) is a 
non- hyperbolic critical point. Most we can say in this 
case is that it is unstable in the directions spanned by 
dy {ey (X exp (3t)), and dz (cz oc exp (6t)). However, the 
phase space drawings (Fig. [2]) reveal that it is in fact 
a unstable critical point, which is the past attractor in 
^p™"'^- the case of a tachyon phantom field living in 
the RS brane, the point (0,0,0) is also the past attrac- 
tor in the phase space vE'^""'^. In both cases this point 
is associated with ricdm = = 0, so that it can be 
identified with an empty (Misner) Universe. The point 
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E = (0, 0, 0) has been included in Tab. HH 

Here we summarize the critical points of the au- 
tonomous systems of ODE ([T6| and (fT7|). and their prop- 
erties (these are also displayed in Tab. ITl| : 

• For the scalar phantom field trapped in the brane 
the critical points are: 

1. The empty Universe solution E — (0,0,0), 
characterized by vanishing matter content 
^cdm = = 0, is the past attractor. This is 
a super-decelerating solution with q — 2. The 
effective EOS parameter of the phantom field 
uJij, is undefined. 

2. The matter-dominated solution M = (0, 0, 1) 
i^cdm = 1, Sl</, = 0) is a saddle equilibrium 
point in 4"^™"^. In this case the expansion is 
decelerated as well {q — 1/2). 

3. The phantom-dominated solution associated 
with the point SP = {-X/V&, ^1 -I- A76) 
{n^ — 1, ujci> < —1), is also a saddle equi- 
librium point in the phase space 5'^™"'^ (the 
third eigenvalue is a positive quantity A3 ~ 
2A2). 



energy-density 




FIG. 4: Energy density vs redshift for the toy model ptot{z) = 
fj,{z + 1)^ -I- (t(1 — exp(— z/(z+ 1)^^)) (Friedmann equation 
(|18l) '). The free parameters have been arbitrarily chosen: 
fi — 0.1, a — 2, 1/ = 0.1. The solid curve represents the 
evolution of the total energy density, while the darker dashed 
line represents the phantom energy density (the cold dark 
matter energy density evolution is represented by the gray 
dashed line). 



• For the tachyon phantom field living in the RS 
brane the equilibrium points are the following: 

1. As in the previous case the empty (Misner) 
Universe E — (0, 0, 0) {flcdm — i^cj, — 0), asso- 
ciated with super-decelerated pace of the ex- 
pansion {q = 2), is the past attractor in the 
phase space vf^T""'^. 

2. The matter-dominated solution M = (0,0,1) 
i^cdm = 1, ^4> = 0) is a saddle critical point, 
characterized by decelerating expansion {q = 
1/2). 

3. The (tachyon) phantom-dominated solution 
TP = (-Ay*/V3,y*,l), where we have de- 
fined = \J {\^ + \/X'^ + 36)/6, is also a sad- 
dle critical point in 4'^''"'=. 

B. Concluding remarks of this section 

With the study of a standard (scalar) phantom model, 
and of a tachyon phantom scenario, we have corrobo- 
rated what we expected. The extra-dimensional brane 
effects modify not only the past asymptotic structure, 
but also the future asymptotic properties of these phan- 
tom models. Actually, by comparing the structure of the 
phase space of the four-dimensional (standard) models in 



the first part of this section, with the one emerging af- 
ter considering the brane effects (former subsection), one 
sees that (see the figures [U [2] and [3]): 



• While there were no past attractors in the phase 
spaces '^ph, *T corresponding to the models ([TU]) . 

respectively, for the same matter content con- 
fined to a RS2 brane, there is a past attractor: the 
empty Misner Universe. Hence, the brane effects 
modify the past asymptotic (early-times) structure 
of the corresponding phase spaces. 

• The matter-dominated solution is not fundamen- 
tally affected by the brane effects. 



• While the phantom-dominated solution was the fu- 
ture attractor in the phase spaces corre- 
sponding to the standard (scalar) phantom model, 
and to the tachyon phantom model, respectively, 
after brane effects are considered, the stability of 
the above solutions is modified: the phantom- 
dominated solution is now a saddle critical point, 
so that, in the corresponding brane models there is 
no future attractor in the phase space. 
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IV. DISCUSSION AND CONCLUSION 

Conventional wisdom convinces us that Randall- 
Sundrum brane scenario, which was proposed as an al- 
ternative to the Kaluza-Klein compactification procedure 
- also to seek for an alternative explanation to the mass 
hierarchy problem -, appreciably modifies gravity only at 
very high energy/short scales (UV modifications only), so 
that it could have impact on the pace of the primordial 
inflation, for instance, but does not modify the late-time 
cosmic dynamics. However, the above statement is true 
only if the energy density of the matter content of the 
brane, dilutes with the course of the cosmic expansion, 
as it is true for most forms of known matter sources (ra- 
diation, dust, quintessential dark energy, etc.). What 
happens if those known forms of matter, are replaced by 
other (exotic) forms which energy density does not di- 
lute, but increases at late times? Take, for instance, a 
Universe filled with a mixture of cold dark matter and of 
a dark energy fluid which energy density grows with the 
cosmic expansion according to ppt = ct(1 — e^""^), so that 
the total energy density of the matter trapped in the RS 
brane is given by: 

The behavior of the total energy density ptot vs the red- 
shift z is shown in the figure ID It is seen that the cold 
dark matter dominates at early times and dilutes as the 
expansion proceeds. Meanwhile, the phantom (dark en- 
ergy) density is a growing function of the redshift and, 
somewhere in the future starts dominating the dynam- 
ics of the expansion. The fate of the cosmic evolution 
(z — > —I) is close to a de Sitter Universe where ptot ~ cr. 
Although this is an unrealistic toy model, several qual- 
itative aspects will be well illustrated. Since the brane 
effects modify the Friedmann equation according to 

then, at very early times where z ^ {2X/ p)^/^ , the Fried- 
mann equation is well approximated by 37?^ = p} jlXc^ . 
This means that the effect of the brane at high energies 
and large redshifts, is to allow for a metamorphosis of 
the background fluid from cold dark matter into a super- 
relativistic stiff fluid. At late times the modification to 
the standard Eisntein's dynamics is also appreciable. Ac- 
tually, at late times, according to standard Friedmann 



equation ^ = ptot ~ c, meanwhile, if brane effects are 
taken into consideration, then, either 37?^ 3cr/2 if a 
is of the same order as the brane tension {a k, A), or 
3i/^ « cr^/2A if cr ^ A. In the former case, aside from 
the factor 3/2, it has to be considered, also, that cr ~ A 
is a very large quantity. In the latter case the energy 
density is boosted by the factor cr/2A which is also very 
large since cr/A ^ 1. This shows the way brane effects 
might affect the cosmic dynamics not only at early times 
but also at late times. 

While the toy model (fT5|) is a unrealistic model with 
no physical motivation at all, other more physically moti- 
vated models, as the scalar phantom, and tachyon phan- 
tom models, also corroborate the former results. Even 
in the case when the effective phantom behavior is the 
result of the peculiar dynamics of a cosmological fluid of 
known nature (the NLED-based model of section II), the 
extra-dimensional brane effects modify the late time dy- 
namics. In all of these models of very different nature, 
the effect of the brane is to modify the stability of the so- 
lution which is associated with the future asymptotics of 
the models: the brane destroys the future attractor, leav- 
ing in its place a unstable saddle equilibrium point. This 
means that the end point of the cosmic evolution is un- 
certain, unlike the case when there is a future attractor 
which is, necessarily, the end point of the cosmological 
evolution. It has to be underlined that the modification 
of the late-time asymptotic structure by the RS brane 
effects occurs only for those exotic forms of matter which 
effective equation of state in GR is less than -1 (phantom 
matter) . 

We want to stress that the above "mixing of scales" 
effect, where both UV and IR modifications of general 
relativity are the result of a same effect (in the present 
case the effect of the extra-space in the Randall-Sundrum 
brane picture), is distinctive only of theories that mod- 
ify the right-hand-side (matter part) of the Friedmann 
equation. In this sense, for instance, it is not to be ex- 
pected that the DGP brane model - designed to modify 
gravity at the IR - can modify GR also at UV ener- 
gies/curvatures. This will be the subject of a separate 
publication llj. 
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